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Abstract 

A study of the order a 2 s corrections to the Drell- Yan (non- singlet) differential cross 
section for incoming states of arbitrary longitudinal helicities is presented. The trans- 
verse momentum distributions, qx, of the lepton pair are studied and the calculations of 
Ellis, Martinelli, and Petronzio (EMP) are extended to include polarized initial states. 
We use the MS scheme and the t'Hooft-Veltman regularization for the helicity pro- 
jectors. From our results one can obtain the bulk of the totally inclusive NNLO cross 
section for the production of a Drell- Yan pair in the non-singlet sector by a simple 
integration over the virtual photon momentum. We show that in the A/5 scheme he- 
licity is not conserved along the quark lines, unless a finite renormalization is done 
and one adapts the physical (MS P ) scheme. This aspect of the calculation is similar 
to the 0(a 2 s ) polarized production of single and double photons. Our spin averaged 
unpolarized differential cross sections agree with the EMP calculations. 



1 E-mail address: schang@phys.ufl.edu 

2 Present Address: Theory Group, Jefferson Lab, Newport News, 23606 VA, USA. 
E-mail address: coriano@jlabs2.jlab.org 
3 E-mail address: rfield@phys.ufl.edu 
4 E-mail address: gordon@hep.anl.gov 
5 Permanent address. 



1 Introduction 



In recent years spin physics has become a very active research area in particle physics. 
This is especially true in Deep Inelastic Scattering (DIS), due to the puzzling results of 
the EMC and SMC collaborations on the nucleon spin and the experimental activity at 
HERA. Various reviews have appeared in the literature discussing the subject in detail and 
highlighting the recent progress made in the area jl], ||, ||| . 

Along with the possibility of probing the spin structure of the hadrons at RHIC and 
HERA, and possibly at other colliders, comes the need to investigate systematically the role 
of subleading effects in polarized collisions. Given the added complexity that spin physics 
brings into the program, it is clear that these new studies will be far more involved and 
will require considerably more effort than in the past, both on the experimental and on 
the theoretical fronts. It is generally expected that radiative corrections to various QCD 
processes in the case of polarized initial states will be substantial and crucial in order 
to reach a better understanding of the behaviour of the polarized parton distributions. 
In particular, in contrast to the DIS case, at hadron colliders the direct gluon coupling 
will allow one to study the size of the longitudinal spin asymmetries of both quarks and 
gluons thereby reaching a more complete understanding of the distribution of the spin of 
the nucleon among its constituents. The Drell-Yan process has received an extraordinary 
amount of attention in the past, both as a way of testing perturbative QCD and as a 
possible mechanism to uncover new gauge interactions at high energy. 

Leading power factorization theorems predict that high energy scattering processes can 
be written in a factorized form, in which a short distance coefficient is convoluted with a non- 
perturbative universal matrix element. At leading twist (or twist 2) these matrix elements 
describe distributions of gluons and quarks in the hadrons, while contributions of higher 
twist describe more complex correlations among partons. In previous works (see for instance 
Q and refs. therein), complete NLO analyses of some important polarized processes have 
been presented. The results show that radiative corrections are sizeable in many cases and 
can be measured in forthcoming experiments at RHIC and possibly HERA. In this work we 
extend the program outlined in by analyzing the behaviour of the Drell-Yan differential 
cross section in the case of longitudinally polarized incoming states. In this work, we present 
a complete NLO study of the non-singlet contribution to the Drell-Yan production of large 
transverse momentum lepton pairs in polarized hadron-hadron collisions, thereby extending 
the unpolarized work of Ellis, Martinelli, and Petronzio Jjj. Our calculations are performed 
in the MS scheme. The factorization and the renormalization at the perturbative order 
are performed in the same scheme. In this scheme, the t'Hooft Veltman Q prescription 
for the regularization of chiral states gives a cross section which is related to the EMP 
result, plus additional helicity non-conserving contributions. In the MS scheme helicity is 
not conserved along the quark lines. By a redefinition of the splitting functions one finds 
that the result for the polarized case can be trivially related to the unpolarized results. 
However, this redefinition results in a different scheme, with polarized splitting functions 
differing from the unpolarized ones by a finite term. This pattern has been already noticed 
in the case of double and single photon production H|. The factorization of the process 
is shown to occur in both schemes, with different structures for the helicity non-conserving 
contributions coming from the real and the virtual diagrams. 
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Figure 1: Drell-Yan process 

The structure of the real and of the virtual corrections are similar to the result of , but 
modified by the helicity dependent factors (1 — h) and h = 4A1A2, where Ai and A2 are the 
helicities of the initial state fermions. This work is mainly concerned with the calculation 
of the radiative corrections to the parton level cross section for the non-singlet case, while 
the hadronic part of the calculation and a phenomenological analysis of our results will 
be presented elsewhere. An extension of our calculation to include the contributions from 
initial state gluons is now in progress. 

This paper is organized as follows. In section 2, we examine the general helicity structure 
of the Drell-Yan process and introduce some of the differential cross sections. In section 
3, we present the calculation of the virtual contributions to the differential cross section, 
while the real contributions are presented in section 4. The study of the factorization of 
the process is done in section 4.1, where we present the final structure of the factorized 
differential cross section. Our conclusions and perspectives are contained in section 5. 

Part of the calculations have been performed using FeynCalc [fic| ], in particular the 
implementation of the Passarino-Veltman [|ll|] recursion procedure for tensor integrals. The 
presentation of the phase space integration can be found in Appendix A, where we describe 
an extension to the polarized case of the previously known methods of integration used in 
the unpolarized case. 

2 The Drell Yan Process 

Deep Inelastic Scattering (DIS) has provided us with a wealth of information on the nucleon 
structure, however, there are limitations on the information obtainable from DIS experi- 
ments in the study of the spin properties of the hadrons. The first and most important 
limitation arises because gluons couple only to higher order. Also, distributions of trans- 
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Figure 2: a s order contributions 

verse spin decouple. For a classification of the various light-cone operators of lowest and 
higher twist see [O], [13] and references therein. The Drell-Yan process is the natural ter- 



ritory to analyze some of the issues connected with these important formal developments. 
For example, it has been observed (see for instance fljl) that the Drell-Yan production of 
charged lepton pairs is a case where transverse and longitudinal asymmetries are compara- 
ble. To leading order, the ratio between the transverse and the longitudinal asymmetries, 
R = Att/All, is estimated to be of order one. This is due to the fact that the gluons 
don't couple to transverse asymmetries and to order a s only the quark-antiquark channel 
contributes. In other processes such as inclusive two-jet production (p p — > j j + X) or 
prompt photon production pp — > 7 + X the gluon channel is available also to leading order 
(in An), thereby making R small. Therefore, it is very important to see how these lowest 
order predictions are affected by radiative corrections. 

Our work is concerned with a study of the differential cross section for the production 
of large transverse momentum lepton pairs, produced by the decay of a virtual photon of 
invariant mass Q 2 . Since the lepton pair is required to have non-zero transverse momentum, 
qx, the differential cross section starts at order a s . At lowest order (LO), one gluon (at 
least) is required to be in the final state. The calculation proceeds up to order a 2 (denoted 
as NLO). 

2.1 Previous Work 

For the unpolarized case previous work on the process can be found in ref. for the non- 
singlet sector. The total cross section, together with a calculation of the K factor and a 
study of the small-x limit has been studied by Van Neerven and collaborators (see [Q and 
refs. therein). In the polarized case, the order a s radiative corrections have been originally 
investigated in [jig] (total cross section). More recently work has been done to extend this 
calculation to the to the Xi?-distribution in |l6j to order a s . The resummation for the of 
the differential cross section at small qr of the lepton pair has been studied in the polarized 
[17]. The integration of our differential cross sections in dimensional regularization 



case m 



provides the bulk of the NNLO (a 2 ) corrections to the total Drell-Yan cross section in 
the non-singlet sector. The remaining contributions, not included in this work, can be 
easily obtained by a simple extension of our work, since they amount to finite (in e) terms. 
However, as we are going to see, the structure of the helicity non-conserving contributions 
do not allow us to conclude that at parton level the total cross section is simply related to 
the unpolarized result by an overall helicity-dependent factor. 
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2.2 General Structure of the Differential Cross Section 

We begin by examining the contribution to the production of virtual photons with invariant 
mass, Q, in hadron-hadron collisions, 

H x + H 2 - 7 * + X, (1) 

from the 2^3 parton subprocess, 

a + b^^ + c + d, (2) 

like those shown in Figs. 4-6. Here H\ and H 2 are incoming hadrons with 4-momenta, Pi 
and P 2 , respectively, and q is the 4- momentum of the virtual photon, 7*, as shown in Fig. 1. 
The 4-momenta of the incoming two parton a and b are labeled by p\ and p 2 respectively, 
and the outgoing 4-momenta of partons c and d are labeled by k 2 and k%. The virtual 
photon 4-momenta is given by fa = q. Conservation of energy and momentum implies that 

Pi +P2 = fa + k 2 + k 3 . (3) 

The hadron-hadron process (|l]) is described in terms of the invariants, 

S=(P l+ P 2 ) 2 , T=(P 1 -q) 2 , U=(P 2 -q) 2 . (4) 

In addition, it is convenient to define the scaled variables 

xt = 2qx/W and r = Q 2 /S, (5) 



where W = \? S is the hadron-hadron center-of-mass energy and qx is the transverse mo- 
mentum of the virtual photon with invariant mass, Q. It is also useful to define the following 
two quantities, 

Q 2 -u lv nr— _ q 2 



xi = - - = & y \Jx T + At, x 2 = ^ — = \e y \jx\ + 4r, (6) 

where y is the rapidity of the virtual photon. The 2 — > 3 parton subprocess (g) is described 
in terms of the invariants, 

s = (pi+ P2) 2 , t = (pi - q) 2 , u = (p 2 - q) 2 , (7) 

and 

S2 = S23 = {fa +k 3 ) 2 , (8) 

where momentum and energy conservation (|3|) implies 

s + t + u = Q 2 + s 2 . (9) 
In QCD the hadronic cross section is related to the parton subprocess according to 

da(Hi + H 2 ^ 1 * + X- W, q T ,y) = 

Gn^aixi, M 2 )dxiGHo^b(x 2 , M 2 )dx 2 I ——(ab — > j*cd; s,t,u)) dtdu, (10) 

V dtdu ) 
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where Gu^aixi, M 2 )dx\ is the number of partons of flavor a with momentum fraction, 
x\ = pi/P\, within hadron H\ at the factorization scale M. Similarly, GH 2 —>b( x 2, M 2 )dx2 is 
the number of partons of flavor b with momentum fraction, x<i = P2/P2, within hadron H2 
at the factorization scale M. In the remainder of this paper we will take the factorization 
scale to be Q and evaluate the parton subprocesses at the same scale Q. Using (||) and (|7|) 
we see that 

s = XlX 2 S, (11) 
(t-Q 2 ) = x l {T-Q 2 ) = -x l x 2 S, (12) 
(u-Q 2 ) = x 2 (U - Q 2 ) = -x 2 x 1 S, (13) 

which from (^) implies that 

X\X2 - X\X2 - X2X1 = T 2 — T, (14) 

where x\ and 22 are defined in @ and 

r 2 = s 2 /S. (15) 

It is now easy to compute the Jacobian 

dx2<it = — — - dx^dy, (16) 

4(xi - xi) 

which when inserted into (|l0|) and integrating over x\ and S2 gives 
da f 1 , f s 2 1 



da f 1 f s 2 

S 2 (W,q T ,y) = / dxi / 
dq^dy J x ™ tn Jo 



where 



and 



dqfdy J x ™ n Jo (xi - x\) 

GH 1 ^a(xi,Q 2 )G H2 -^b(x2,Q 2 ) s^—(ab^7*cd;s,t,u), (17) 

dtdu 



xix 2 + t 2 -t 

x 2 = z , (18) 

X\ — X\ 

A = (t-x 1 + x 1 (1-x 2 ))S, x ™™= d ^L. (19) 

1 — X2 



The maximum value of r 2 arises when X2 = 1 in ([14]), while the minimum value of x± occurs 
when T2 = and X2 = 1. 

For 2 — ► 2 parton subprocesses such as the Born contribution in Fig. 2 one has 

s J?-( s ' M ) = ^2) s-?(o + 6 ->■ 7* + c; s, t), (20) 
atcm at 

which when inserted into (17) results in 

da f 1 , 1 



s i~2~r( w ^T,y) = [ dxi 

dq&dy J x ™ m 



dqj.dy J x ™ in (x\ - x{) 

G Hl ^ a (x 1: Q 2 )G H2 ^ b (x 2 , Q 2 ) s^(ab -» 7 *c; s, i), (21) 
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where in this case s 2 = and s + t+u = Q 2 . Finally, the Drell-Yan differential cross section 
for producing muon pair with transverse momentum q^ and invariant mass, Q, at rapidity 
y in the hadron-hadron collision, 

Hi + H 2 - (7* - M + /0 + X, (22) 

at center-of-mass energy, W, is given by 

' " A S- J ^ r (if 1 ff 2 ^7* + A-;H',, 7 ., !/ ), (23) 



\3irQ 2 J dq^dy 

where a is the electromagnetic fine structure constant and where the virtual photon dif- 



ferential cross section is given by ( |17| ) or (21). In the remainder of this paper we will 



concentrate on the virtual photon differential cross section. 
2.3 The Non-Singlet Cross Section 

We define the non-singlet to be the difference of the antihadron-hadron and the hadron- 
hadron differential cross sections as follows: 

dqj,dy 

S ^r( H + H 2 ^ 1 *+X)- S-^—(H + H 2 ^ 1 * + X). (24) 
dqj,dy dqj,dy 

It is easy to see that 

d<7NS rs™ ax ^ 

dqfdy fi, Jo [Xi - 



f,f 

iNS /„ ^2\r-iNS ^>2\ _d&NS , 



Xi) 



G%t qf (xi,Q 2 )G%t qf ,(^Q 2 ) s ^( s ^ u )' ( 25 ) 
where the non-singlet parton-parton differential cross section is given by 

and the non-singlet structure function is defined according to 

G^^q = Gh->q — Gn^q- (27) 
The sum in (25) is over the quark flavors / and /'. As we will see this sum turns out to be 



diagonal in flavor space so that one needs only to sum over / = /'. 

We will organize the Feynman diagrams in the same way as was done in J?J. At order 
a s the Born term comes only from the square of the two amplitudes in Fig. 2, which we 
write schematically as 

da B ~ \Li + L 2 \ 2 . (28) 
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The Born term contributes only to the diagonal part of the non-singlet differential 

cross section in (^5|) since it is non zero only when / = /'. 

The order a 2 , contributions are separated into the real and the virtual diagrams. The 
virtual term, do mrtual comes from the interference between the amplitudes L\ and L 2 and 
the eleven amplitudes shown in Fig. 6, 



do 



virtual 



2Re^(L 1 + L 2 )J2vA , (29) 



which again contributes only to the diagonal qjqj part of the non-singlet differential cross 
section in (25). 



The other real diagonal q/qj diagrams are 



i=l i=\ 



where the amplitudes Fi are shown in Fig. 3 and the amplitudes Gi are shown in Fig. 4 
with / = /'. The absolute square of the Fi amplitudes can be written as follows: 



8 

|2 



i=i 

I 4 8 \ 8 

l^i + F 2 \ 2 + \F 3 + F 4 | 2 + 2Re ]T F £ F* + | £ F\ 2 , 

\i=l i=5 / i=5 



(31) 



where we have used the fact that 

Re{{F 1 + F 2 )(F 3 + F 4 )*) = 0. (32) 
Thus the diagonal real q~fqf contributions are, 

d ° T q e f a q f = d ° l + d ° 2 + d<5 "3 + da F , (33) 

where 



den ~ \F 1 + F 2 \ 2 + \Y / G i \ 2 , (34) 

i=l 

do 2 ~ l^s + ^l 2 , (35) 
da 3 ~ 2Re(j2F i Y / F*\, 

\i=l j=5 / 



(36) 



d<j F ~ lE^i 2 - ( 37 ) 

The real diagonal qjqj diagrams are 

j=l j=4 
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where the amplitudes Hi are shown in Fig. 5 and the relative minus sign between the direct 
and exchange diagrams is due to Fermi statistics since / = /'. The absolute square of the 
Hi amplitudes can be written as follows: 



i=l j=4 
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4 8 / 4 8 \ 

i=l i=4 \i=l i=4 / 



(39) 



Thus the diagonal real (//<?/ contributions are, 



if if 
where 

l 



da r Q e f ^= da H + da A , (40) 



da H ~ |^F,| 2 + |^^| 2 , (41) 

i=l i=4 

d(7 4 ~ -2ifefc (42) 

\i=l i=4 / 



The diagonal (/ = /') real part of the non-singlet cross section in ( pq ) is given by subtracting 
(pC|) from (|3^) as follows: 

da^s = d&i + da 2 + da 3 - da A , (43) 
where we have used the fact that d&F = dan which arises due to 

Jps 3 ^ Jps 3 fr[ Jps 3 fri 

When integrating over the two identical particles k^ and k A over the phase space PS% in 
(p-41|) an extra statistical factor of 1/2 must be inserted so as not to double count. Thus, the 
diagonal (/ = /') part of the complete non-singlet cross section in ( |26| ) (real plus virtual) 
is given by 

d&Ns = da i + d&2 + das — da A , (45) 

where 

da x = da( act + da virtual , (46) 
and where da( act is the "factorized" cross section given by 

da{ act = da t + da1 ounter . (47) 

It is necessary to add a "counterterm" cross section ( j i ^ ounter to da± to cancel the initial 
state collinear singularities. 
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The off-diagonal real qfqp diagrams are 



4 



^ -lE^il 2 . ( 48 ) 



i=i 



where the amplitudes Fi are shown in Fig. 3 and where / ^ f. Similarly, the off-diagonal 
real diagonal qjqf diagrams are 

i=l 

where the amplitudes Hi are shown in Fig. 5. In this case, 

|£f,| 2 = iI>I 2 , (50) 

i=l i=l 

so that these two contributions cancel when forming the non-singlet cross section in (p6|). 
Thus to order a 2 (^) becomes 

dPNS ftitr \ ST^ I j / 2 j ^ 



s j 2~r ( w iQT,y) = J2 dxi ds 

dq^dy J x ™n Jo 



(xi - Xi) 



/ 

G N H l q (x^Q 2 )G N H s (x 2 ,Q 2 ) s^(s,t,u), CD 



nf ^ L ,^ ^H^qjm-* > J dfdu 

where the sum is now only over the quark flavor / and where 



d& N s da B , dai da 2 da 3 da 4 

s = s-r— — h s—— — h s— - — hs-ri s , , . (52) 

dtdu atdu atau dtau atau dtau 

2.4 The Helicity Dependence of the Born Term 

The lowest order (non-singlet) contributions to the large transverse momentum production 
of virtual photons (see Fig. 2) with invariant mass q 2 = Q 2 arise from the two, q+q — > 7* +g, 
"Born" amplitudes L\ and L 2 - We refer to these two diagrams as the direct and the exchange 
(or crossed) amplitudes, respectively. The absolute square of the sum of these amplitudes 
form the Born differential cross section in (|28|). In Fig. 2 we have also generically illustrated 
the expansion of the amplitude which appears in da/dQ 2 up to order a s . Notice that the 
quark form factor contributions in Fig. 2 (and the related 0(a 2 ) corrections, not included 
in the picture) do not appear in the study of the cross section, da /dq^dy, for transverse 
momentum, qT, greater than zero. The two-loop 0(a 2 ) corrections to the quark form factor 
can be added in the study of the total cross section, da/dQ 2 , by interfering the 2-loop on- 
shell quark form factor of ref. [18| with the lowest order qq'y* annihilation channel and by 



using helicity projectors for the initial quark states. 

We calculate the spin dependence of the cross section by using the helicity projectors, 
P± = ^(±±75), which project out the helicity states of an initial state quark and antiquark, 
respectively, as follows: 

1 1 _ 

u(pi,/ti) = -(1 + hi^ 5 )u(pi), v{p 2 ,h 2 ) = -v(p 2 )(l + /1275), (53) 
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where hi = ±1 = 2Ai corresponds to quark helicity Ai = ±|, and h 2 = ±1 = 2A2 
corresponds to antiquark helicity A2 = 

The squares of the direct amplitude L\ and exchange amplitude L 2 in Fig. 2 in iV = 4— 2e 
dimensions are given by 



Mdd(h\, h 2 ) = e} 9 2 9 2 s ^ ((1 - ef - h^l + e) 2 ) , 
M cc (hi, h 2 ) = e 2 f g 2 g 2 ^f- ((1 - e) 2 - h x h 2 (l + e) 2 ) . 



(54) 



where a = <? 2 /47r is the fine structure constant, and e/ the charge of the quark. The 
quantity Cp/N c is the color factor, and a s = g 2 /4:ir is the QCD strong coupling constant. 
We pause here to note that in our NLO calculation a s — > a s (n 2 ) is the running coupling 
constant renormalized at the scale fi 2 in the MS scheme satisfying the renormalization 
group equation 



2 da s 



— OLs 



a s 



47T 



a s \' 



47T / 



(55) 



with 



A) 



11N C 2Np 
3 3~~ 
34iVg iqn c n f 



2C F N F 



(56) 



and with Nf the number of active quark flavors. The interference term is more complicated 
and is given by, 



2M dc (h 1 ,h 2 



2 2 2 Cf 4 

e f 9 9 °W r Tu 



(1 - e)(Q 2 s - etu) -hih 2 (l + e)(Q 2 s + etu) - 2h 1 h 2 etu 

The sum of the two Born amplitudes squared is 

\M B (h u h 2 )\ 2 = \M B {h)\ 2 = M dd (h) + 2M dc (h) + M cc (h) 

= *)?&%Yu [ (1 - e) ( 2Q2s + (1 - e){t2 + u2) ~ H 

-h(l + e) (2Q 2 s + (1 + e){t 2 + u 2 ) + 2etuj - Ahetu , 

and depends only on the product h = h±h 2 . 

To any order in perturbation theory we can write, 

|M(/n, /i 2 )| 2 = \M(h)\ 2 = (1 - /i)|M| 2 + /i|M ++ | 2 = |M| 2 + hA\M\ 2 , 

where h = h\h 2 = 4A1A2 and where 

\M\ 2 = \\J2 M{h u h 2 )\ 2 = W \M(hi,h 2 )\ 2 , 



(57) 



(58) 



(59) 



(60) 



hi,fi2 
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is the spin averaged (unpolarized) amplitude squared. Furthermore, 



|M_ 



M ++ | 2 and |M_+| 2 



I M 4 



so that 



\M\ 2 = l -(\M ++ \ 2 + \M + 4 2 ). 
The spin asymmetry, A|M| 2 , is defined according to 

A|M| 2 = \ (|M ++ | 2 - |M + _| 2 ) = |M ++ | 2 - |M| 2 . 



(61) 
(62) 

(63) 



The spin averaged (unpolarized) amplitude squared is determined from \M(h)\ 2 by set- 
ting h = with A|M| 2 the coefficient of h. For the Born term, this results in 



i I eW^l [ (1 " £) ( 2Q2s + (1 " e){t2 + ~ 2etU ) 



and 



A|M| 2 : 

Adding these two terms yields 



222 Cf 2 

e f g 9 °W c Tu 



(1 + e) (2Q 2 s + (1 + e)(t 2 + u 2 ) + 2etu) + Aetu 



\M. 



++I 



2 2 2 C.F 8e 

" e / 5 9 °W c Tu 



Q 2 s + (t + uf 



(64) 



(65) 



(66) 



which is proportional to e and vanishes in the limit e — > 0. Since at the Born level there 
are no ^ singularities that might combine with this term to yield a finite contribution, in 
the limit e — > 0, 

|M ++ | 2 = and A|M| 2 = -|M| 2 . (67) 

For the q + (/ — > 7* + g subprocess the condition that |M ++ | 2 = means that the quark 
helicity is maintained (does not flip) in the collision. The incoming quark line with helicity 
±i turns around and becomes an incoming antiquark with helicity =F |- We refer to this 
as "helicity conservation". We see that at the Born level helicity is conserved in the limit 
e -► 0. 

In N = 4 — 2e dimensions the differential cross section is related to the 2^2 invariant 
amplitude according to 



da 1 /MttsV 1_ 

S dt (S,t ' /l) " 16^ \tu~) r(l-e) 



\M(h)\ 2 , 



which can be written as 



or as 



da , , . . , . da^ , . , da 11 , 
s— (s,t,h) = (l-h) s-jf(s,t) + h s-±±(s,t), 



(68) 

(69) 
(70) 
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where sda^/dt is the unpolarized cross section 



dcP 1 / da 



dt 



dt 



+ s 



da + - 
dt 



and 



At the Born level we have, 



da LL _ 1 / dd++ 
dt "2 S (If 



— s- 



da^ 

IT 



s ^ {s ,t) = e}K 2 ^T B (Q*,u,t), 
dt 1 s 

s^j^(s,t) = -e 2 f K 2 ^-eA B (Q 2 ,u,t), 



where 



T B (Q 2 , u,t) = — [(1 - e) {2Q 2 s + (1 - e)(t 2 + u 2 ) - 2etu 
= 2(l-e) 



t u 



ut 



and 



A B (Q 2 ,u,t) = — Q 2 s + (t + u) 2 
tu L 



where K 2 is defined by 



Ko = Tta 



C h 



1 



4 ¥ 2 sQ 2 



N c r(l -e)\ Q' 



tu 



where we have rescaled, a s — > a s (^ 2 ) e , so that it remains dimensionless in 
dimensions. At this order of perturbation theory, we have 

s^-(s, t, h) = e 2 K 2 ^ [(1 - h)T B (Q 2 , u, t) - heA B (Q 2 , u, t) 
As we saw earlier, in the limit e — ► helicity is conserved so that 

s^(s,t,h) = (l-h)s^-(s,t), 



which implies that 



t — —(s,t) = and s — ?-(s,t) = -s—^-(s,t), 
dt y ' dt 1 ; dt y h 



so that to leading order, 



2 da LL Ida 11 

ALL = ^rhr = ~ l - 
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To connect with the notation of ref . |?J , we note that 

K 2 T B (Q 2 , Ul t) = KT (Q 2 ,u,t), 
where K and To are define in ref. @ as 

C F (1-e) / 4ith 2 Y ( sQ r 



K = 2vra 



N c r(i - e) V Q 2 



tu 



and 



T (Q 2 ,u,t) 



t uj ut 



(81) 



(82) 



(83) 



2.5 General Structure of the Helicity Dependent Cross Section 

To order a 2 the non-singlet helicity dependent Drell Yan cross section can be written 



according to (|51|). Namely. 



ag^ay V ix™» Jo (xi - 



Xl) 



da NS 



dtdu 



(84) 



where Ai and A2 are the helicities of the initial two hadrons, respectively, and h = 4A1A2 
and where the sum is now over the quark flavor / and the parton helicities Ai and A2. The 
non-singlet polarized structure functions are given by 



H[A]^q f [X] ~ ^H[A]^q f [\] ~ ^ H [K]-*q } [A] ■ 

The parton-parton non-singlet differential can be written as a function of h as follows: 

daws , 



(85) 



dtdu 



(s, t, u, h) 



da 4 
s - - (h). 
dtdu ' 



(86) 



Prom ( |7q ) we see helicity is conserved for the Born term so that 

s^j~(s,t,u,h) = (1 - h)5(s 2 ) s^-(s,t), 
dtdu dt 

where da^ is the unpolarized cross section in (|73[). Helicity is also conserved for da 2 defined 



(87) 



in (|3^) and for d&% defined in ( p6| ) so that 

da 2 



da% 



and 



(s,t,u,h) = (1 - h) s— — (s,t,u), 
dtdu dtdu 



s——(s,t,u,h) = (l-h) s——[s,t,u). 
dtdu dtdu 
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The Born term and these two cross sections have the property that dfrr + +) = as 
expected from helicity conservation. The term d&4 defined in (|42] ) is different. Here helicity 
conservation implies that da^ H — ) = and we have 

11(74 (s,t,u,h) = (l + h) s^(s,t,u). (90) 



dtdu dtdu 

The term da\ is complicated and whether helicity is conserved or not depends on how one 
does the factorization in (47). We write the helicity structure as follows: 

ddi , , s , , x daV , , , ddn at , , , 

W S '*'^ = (1 - h) s d± {s ^ u) + h s *k is ^ u) > (91) 

where daQ at is a regularization scheme dependent helicity non-conserving piece and daf is 
the unpolarized spin-averaged cross section. 

To order o? s the non-singlet parton level spin asymmetry defined in (72) becomes 

da{h , . da^Q dan at daf , . 

s^^(s,t,u) = -s—^ + s—^- 2s— (92) 

dtdu dtdu dtdu dtdu 

where da^ s is the non-singlet unpolarized spin averaged cross section given by 

da^a , . da^ daf da? da? daf , , 

atau atau atau atan atau atan 

At the parton level it is no longer true that All = — 1 as was the case for the Born term 
in (180). 



3 Virtual Diagram Contributions 

The list of the diagrams with virtual corrections contributing to the non-singlet sector is 



defined in (2£) and shown in Fig. 6. We have omitted all the self-energy insertions of quarks 



and gluons and the ghost contributions. 

Our calculations are performed in the MS scheme using dimensional regularization 
to regulate both the ultraviolet and the infrared singularities. We remove the ultraviolet 
singularities in the relevant subdiagrams by off shell regularization. Then by sending on 
shell the initial state quarks and the final state gluon, we encounter singularities in the 
form of double poles and single poles in e = 2 — n/2. The renormalization of the diagrams 
is enforced in such a way to get straightforwardly final results free of UV singularities by 
symbolic manipulation. We briefly comment here on how this is achieved in our case, 
although a more detailed discussion of the method will be presented elsewhere. As a 
first step, we start with the Passarino-Veltman reduction of the tensor contributions and 
generate the expression of the invariant amplitudes of the reduction. In our case only one 
external line (the photon ) is massive. The calculation is ambiguous because of the presence 
of massless partons and, in principle, one has to be particularly careful to differentiate 
between the UV and the IR origin of all the singularities, which in our approach happen to 
be identified {ejjv = e iR i n Dimensional Regularization). Therefore, the expressions of the 
invariant amplitudes in the Passarino Veltman reduction are, at this stage, unrenormalized. 
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At the second step, we renormalize the scalar diagrams which appear in the result by 
removing the UV poles, (since these are removed by coupling constant renormalization) 
and then switch e — > — e in the final result. This allows us to at this stage interpret the 
left over singularities as IR singularities. Even at this stage the procedure is not complete. 
One needs a prescription to handle the massless tadpoles (self energies at zero momentum) 
in a consistent way. The Passarino Veltman procedure in the massless case, in fact, suffers 
from ambiguities of this type, since massless tadpoles are intrinsically ill-defined (being 
scaleless). In dimensional regularization they are usually set to vanish. However, since we 
have identified UV and IR singularities in our algorithm, this last step is also ambiguous. 

At the third step we impose two prescriptions which allow us to handle consistently way 
all the tadpoles and all the massive self energy contributions. These include : 

1. isolated tadpoles (B(0)) and "linear" tadpoles (or eB(0) terms). 

2. massive self energy contributions (B(q 2 ) and eB(q 2 )). 

In (2) we implement off-shell renormalization, while in (1) we reinterpret all the scale- 
less contributions as the massless limit of the off-shell ones. After step 3, we obtain the 
correct renormalized expressions of all the Passarino Veltman coefficients for all the one- 
loop amplitude. The method works at one loop for any n-point function and is extremely 
convenient in practical applications since it allows us to forget about the nature of the 
singularities after IR/UV identification. The symbolic implementation of the method is 
also quite straightforward and will be illustrated elsewhere pOfl . 

The one- loop order result for the cross section for the process q + q — > 7* + g is given by 

d<r virtual , 

-(s, t, u, h) = 



dtdu 

e}K 2 ^-5(s + t + u - Q 2 ) I ((1 - h)T B - heA B ) 



a s T(l 



2vrr(l -2e) 



2 



+ — ( 3C p — 2Cp In — ^ + —Nc + Nc In 
2 , , 2Q 2 s + t 2 + u 2 Q' 2 (t 2 + u 




tt z (AC f + N c )— — ■ 2(2C F - N C )- 



3tu tu(t + u) 

8(2Q 2 s + t 2 + u 2 ) Q A s(t + u) t 2 + u 2 \ 

t) (Q2- u )(Q2-t)) 

N c 2 -^ + 2c/ + {s + U ^ 
tu J 

+ 2c/ + is + tr 




+2{2C F - N c ) \ Li 2 

Q 4 - (i + u) 



t tu 
t + u\ 2Q 2 s + u 2 + t 2 + 2s 2 




(t + u) 
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(94) 



where h = h\h2 = 4A1A2. This result agrees with the (unpolarized) spin-averaged (h — > 0) 
result given in |7]|, although its structure is more simplified since we have performed some 
additional manipulations on the finite contributions. Therefore, we confirm the result of 
Q, which had been used before by other authors 21, 22], but never checked before. The 
quantities Tg and As appear in at the Born level and are given by equations (|74|) and 
(frq) , respectively. Notice that the double poles and single pole structure (in the h —* 
limit) automatically reproduces the singularities of the virtual contributions of ref. In 
particular, the double poles in e correctly multiply the two-to-two Born contribution Tg, 
which is the Born level unpolarized result. The presence of the heAs term in (93) implies 
that in the t'Hooft-Veltman H), scheme the virtual corrections by themselves do not 
conserve helicity (i.e. they do not satisfy the equation (|78|)). In particular, the finite part 
of da++ 1 dt in the limit e — > for the virtual corrections is given by 



^virtual 



ejK 2 ^-A B (s,t)[3C F 



'2vrs 



11 

IT 



2CWn-2 +— N c + N c ln 



ut 




(95) 



where s + t + u = Q 2 . In some regularization schemes, such as those enforcing an anti- 
commuting 75 in n dimensions, helicity would be conserved at parton level. In the t'Hooft- 
Veltman pi 19] scheme helicity is not manifestly conserved. 



4 The Real Emission Processes 

From (|5l|) and (52) we only need to consider real contributions that are diagonal in flavor 



space. The diagonal (/ = /') real part of the non-singlet cross section in (26) is given in 
© by 

doxs = da 1 + da 2 + da 3 - d& 4 , (96) 



where 



d&! ~ \F 1 +F 2 \ 2 + \J2 G 



2 

1 1 > 



(97) 



i=l 
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and where d&2, da^, and da 4, are defined in p5|), (p6|), and (^), respectively. The cross 
section da\ requires factorization, 

da{ act = da l + da1 ounter , (98) 

where the "counterterm" , dfr<x> unter to da±, cancels the initial state collinear singularities. 
The factored cross section is then added to the virtual term to give, 

da x = da{ act + da viTtual , (99) 

as in (f26|). 

The reason why we include the two diagrams Fi and Fi in this partial contribution 
is because it is only after adding these two contributions to the Gj diagrams that the 
corresponding cross section factorizes. Notice that the two diagrams G4 and Gs generate 
collinear singularities which are different from those coming from the remaining diagrams 
of the same set. In fact, while in most of the G-diagrams the gluon can become collinear to 
an incoming quark line, in G4 and Gs we have a gluon splitting into two gluons which can 
become collinear. The latter singularity is characterized by the invariant S23 — ► 0, and it 
appears also in the two diagrams denoted F\ and F2 from the F-set. These two singularities 
do not require factorization since they are directly canceled by corresponding ones in the 
virtual contribution. The virtual diagrams responsible of this cancelation involve the gluon 
self energy corrections to L\ and L2 in the final state, and the emissions of gluons, ghosts, 
qq and q'q' pairs in the final state, as shown in Fig. 7. As noted by EMP, all the explicit 
singularities in the real emission processes come from this mixed F-G set of diagrams. 
We use the word explicit because there are other singularities due to collinear emission 
in the remaining sets F5,...,F% and Hi,...,H$ but they cancel out when the non-singlet 
combination is formed and need not be f actor ized. 

In a scheme with a totally anticommuting 75 we would clearly get a selection rule for 
the helicity contributions of the form, M ++ = 0. In fact all the final state emissions are 
vector-like and an incoming antiquark of positive helicity can be replaced by a quark of 
negative helicity flowing back toward the initial state. 

If it is postulated that 75 anticommute with the other 4-dimensional Dirac matrices and 
anticommute with the remaining ones, then all the corresponding algebraic relations can 
be shown to be consistent with dimensional regularization Therefore we define 

757 M + 7 M 75 = 
757m - %75 = 0. 

(100) 

Application of this definition of 75 leads to the subdivision of all the kinematic momenta 
into 2 sets : 4 dimensional ones and n — 4 dimensional ones, the latter called hat-momenta. 
All the dependence on the hat-momenta can be removed from the virtual corrections. 
The hat-momenta are removed only after performing the internal loop integrations using 
dimensional regularization. Notice that this procedure introduces spurious finite terms 
which survive even in the limit in which e — ► (n = 4 — 2e). These terms are due to 
the fact that single infrared poles in 1/e can be canceled by e-dependent factors generated 



18 



by the trace in the n — 4-dimensional subspace. There are other special features of this 
regularization which deserve special attention. 

In the MS scheme, the HVBM prescription for the regularization of the chiral states 



generates the helicity non-conserving term shown in (91). In general one would expect 
that all the infrared sensitive contributions in the real emission diagrams, and therefore 
the G, the F and the H diagrams, would similarly contribute to an helicity non-conserving 
term. However, as we are going to show in the next section, this is not the case. The issue 
of the helicity conservation is completely settled by the explicit singularities and by the 
diagrams of the first set {da\). This result suggests that the helicity non-conservation is not 
connected to all the singularities, but only to those singularities which require factorization. 
This observation is crucial in order to simplify the final result for both the virtual and the 
real contributions and identify their factorized form. The factorized do\ cross section can 
be written in the form 

da{ aCt , ,n dam, , dak at . . 

S ^7 (s ' *' u ' h) = s W s ' *' u ' h) + h s W s ' *' n) ' (101) 

and combining with the virtual contribution gives 

dd\ . . dd\R , . , d<7n at ,* nn s 

s-r—[s, t, u, h) = s——{s, t, u,h) + h s—f—{s, t, u), (102) 
atau atau atau 

where 

d&m, ,n da 1R . daf rtual . 

S Hd^^ *' U ' k) = S lM^^ *' U ' k) + S ^T (S ' *' U ' h) - (1 ° 3) 

In the HVBM prescription both d(T\R and da\ irtual are helicity non-conserving, but the sum 
dam conserves helicity. Namely, 

where daf R is the spin-averaged (unpolarized) result. Notice that since da hat is proportional 
to h it does not contribute to the spin-averaged cross section which comes from setting 
h = 0. 

We give the result for the non-singlet cross section in two different schemes: the MS 
scheme and the so called MS P scheme Q. In this later scheme we subtract additional 0(e) 
contributions to the polarized AP qq splitting function, resulting in 

AP qq (z) = P qq (z) + 4e(l - z). (105) 

In this scheme helicity conservation at the parton level is maintained (d&Q at = 0) but 
P qq AP g q. In the MS scheme P qq = AP qq and helicity is not conserved at the parton 
level (daQ at 7^ ). Since we move from one scheme to the other by a finite renormalization, 
the difference between the two evolutions is absorbed in the 0(a s ) (i.e. Z^) contribution 
to the renormalization factor in front of the parton distributions. These changes imply 
that the NLO contribution to the kernel (i.e. P^) is also modified by an amount — > 
P^ — 0qZ^\ where (3q is the first coefficient of the beta function (see sect. 10.2 of ref. [23| 
for more details). 
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4.1 Factorization of the Real Contributions and da\ 



The standard procedure to be used in the calculation of the NLO contributions is to add the 
real emission to the virtual contributions, thereby canceling those singularities which are 
characterized by double pole in 1/e. After factorization of the collinear singularities in the 
initial and final states the final result of the calculation is expressed in the form of regular 
terms and of plus distribution in the the variable s 2 in (||). The two "plus-distributions" 
in question are 1/{s2)a+ an d (}n(s2) / s 2 ) a+ , which are defined in terms of the upper limit 
of integration in (|l7|) and defined in (19). These two distributions are defined by 



ds 2 f(s 2 ) 



(S2)A+ 



A 1 



lns2 

S2 



S 2 

A lns2 



(/( S2 )-/(0)), 







S2 



(/(-a) - /(0)) • 



(106) 



As usual we define the relation between the bare and the renormalized structure functions 
by the equation 



G A ^(x,M 2 



- 1 dz_ 

z 

We choose for R the usual form 



&ij5(z 



Gbare 



(107) 



R 



e JK J T[l - 2e] 



M 2 



(108) 



where the finite pieces C^Az) are arbitrary. In the MS scheme they are set to be zero. 
The polarized splitting functions, APy, are given by 



AP, 



AP 



C F 
N F 
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2 

AP 99 = N C 
AP gq = C F [2- 
with the distribution 1/(1 



(l + z 2 
(l-z)_ 

[2« - 1] 



1 1 



defined by 



(1 



(1 



1 



+ 



33 - 2N F 
6 



5(1 



1 /(«)-/(!) 



(109) 



The factorized cross section is given to order ct^ by the renormalized cross section da± 
with the the collinear initial and final state singularities subtracted 



- fact 



dtdu 



da\ j 
dtdu 
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r 1 da 
J2 J o dztRk^zuM^s-^- \ pl ^ ZlPl 6[ Zl (s + t-Q 2 ) + u) 

5(z 2 (s + u-Q 2 )+t). (110) 

^ Jo at 



The collinear and soft divergences of the two unobserved jets are described by the limit 
S2 — ► and will appear dimensionally regulated in the form 1_<E . These divergences can 
be expanded in terms of poles in 1/e and regulated distributions by the identity 



1 1 x( \ 

5(s 2 ) 



1 - elnA + -e 2 In 2 A 
2 



1 /lns 2 \ 
+ — e (111) 

\ S2 /, 



(S2)A+ \ S 2 J A+ 



up to order e 2 . The collinear singularities associated with the unobserved jets cancel 
when the contributions from the three quark-antiquark annihilation processes qq —* gg^f* , 
qq — ► q'q'^f* and — ► 997* are added together. The soft component cancels against the 
corresponding poles from virtual contribution. Adding real and virtual contributions and 
factorizing the remaining collinear singularities with eq. ( |110| ), all the divergences cancel 
and the result is free of all singularities. 

Details of the three particle phase space integration are given in appendix A. In all cases 
the polarized and unpolarized three-body matrix elements are integrated simultaneously 
since, except for the presence of the hat-momenta present for the set of diagrams labeled 
G and F\ and F%, they have a similar structure. 

The diagrams | J2t C«| 2 an d l-^i + -F2I 2 are integrated in 4 — 2e dimensions since they 
generate soft and collinear, and collinear only poles respectively. The non-4-dimensional 
parts of these matrix elements are different in the polarized and unpolarized cases in the 
HVBM scheme. This leads to a different structure for the integrated three-body matrix 
elements. In both cases there are also finite contributions from the hat-momentum integrals. 
In the case of the diagrams |iq + -F2I 2 the polarized and unpolarized results differ by 

x akt u ,Q 2 s + 2t 2 + 2tu + 2u 2 

5a ~ 4Nfd(s2)— (112) 

where external color factors and coupling constants have been neglected. This difference 
comes both from the 0(e) part of the three body matrix element when multiplied by the 
factor — 5(s 2 )/e coming from the expansion of s^" ~ e as well as from the hat-momenta. 

Similarly for the | J2iGi\ 2 terms there are differences between the polarized and unpo- 
larized result arising from the differences in the matrix elements up to 0(e 2 ) as well as the 
hat-momentum contributions. In this case the differences coming from the O(e) part of the 
matrix element are too complex to present here. The total contribution coming from the 
hat-momenta in the MS is given by 

da hat , , dak at , . , . 

S ^ {s ^ h)=hs ^ {s ^ (U3) 

where 

dvn at , \ 2 T , ol 2 \ 11 Ar , (t + u) 2 n „ 
dtdu TTS o tu 



u-3Q 2 



u(Q 2 S2 - tu) 
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J 1 ( 2s 2 t + Q 2 (s 2 + t) 2Q 2 ( 



+ 



(S2 - ty \ 
i 



Q 2 s 2 - tu 



t 2 



Q 4 + u(2s 2 + u) - 



2s 2 t 



U \ (s 2 ~u) 2 , 

2Q 2 s 2 (2uss 2 + Q 2 (t-u) 2 ) 

u(s 2 - U) 2 



+ {t <-> u} 



[1U) 



When the three-body contributions from \Fi + F 2 \ 2 and the | J2 i G{\ 2 diagrams are all 
added together and the collinear poles factorized according to Q11C| ) then the result is given 
by 



s^(s, t, u, h) = e)K 2 ^ J ((1 - h)T B - heA B ) 6(s 2 ) 
dtdu J s I 

(2C F + N C 1 . s 11 



as r(l - e) / 47r/i 2 
27rr(l-2e) I Q 2 



+ — I 3C F — 2C F In 



~-, + —N c + N c ln 
Q z 6 



ut 



-N F 



+- 



Ots 

'2^ 



11 



N c 



1 



Q 2 



3^ ln T + 



67 



7T 



3 +ln (^) 



+ (2C,-iV c )|2f^^) 



(s2)a+ 



V «2 / A- 

2, 



+ 



1 



ln 



~2N F - 11N C 

. 6(s 2 )a+ 



+e 2 i4: 2 ^(l-h) 



ut - Q z s 2 
(u - s 2 )(t - s 2 ) 

N c 
2 



+ 4C F 



(s 2 )a+ \(u - s 2 )(t - s 2 )(ut - Q 2 s 2 ) _ 
ln( S2 /Q 2 )^J "j 



s 2 



11 s + Q 2 s 2 (3s 2 -4t) 



2s 



+ 



Q 2 



6 t-u 2tu(t-s 2 ) 2 u(t-s 2 ) 3t 2 



c | - g | 2Q 2 n + ts 2 -2Q 2 s 2 | 2Q 2 (t-u) 

\ (t — s 2 ) 2 t — s 2 ut t(tu — Q 2 s 2 ) tu(t — s 2 ) 



2\2 



+C F ln 



Q 2 \ Aut(u - Q 2 ) - 4Q 2 (u - Q 2 ) 



s 2 



ut(ut — Q 2 s 2 ) 



uts 2 Q 2 — u 



(t ~ S2) 2 



+ 



2Q 2 



t 2 Q 2 \ . 
t(t -s 2 ) t t 2 K ' 



+N F 



2stu -Q 2 (t- u) 

sl 2 ^ 2 



+ 2(2C F - N c )hx 



ss 2 



Q 2 + s 



+(2C F -N c )ln 



ss 2 



{u - s 2 )(t - s 2 ) J la 
(Q 2 -s)((t-Q 2 ) 2 + (u-Q 2 ) 2 ) 



^on i ( ut ~ Q 2g2 ^ 
+2C ^[(u-s 2 )(t-s 2 )) 



(ut - Q 2 s 2 ) ) (t — s 2 )tu(u - s 2 ) 

Q 4 + s 2 + (Q 2 -t) 2 + (Q 2 -u) 2 



(tu — Q 2 s 2 )tu 



where 



da 



~ fact 



dtdu 



(s, t, u, h) 



dam . , . dak at . . 

s - , (s, t, u, h) + s (s, t, u, a). 



dtdu 



dtdu 



(115) 
(116) 
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Notice that the term that survives when h = in da\R exactly cancels the corresponding 
term in davirtual, resulting in the helicity conservation implied by ( |104| ). 



4.2 The Finite Cross Sections dcr 2 , da 3 , and da 4 
The remaining real contributions to the non-singlet cross section, 

da 2 ~ \F 3 + FA 2 , 



da 3 ~ 2Re(Y,FiJ2 F i) 

\i=l i=5 / 

8 

ia F ~ lE F *| 2 ' 



(117) 
(118) 

(119) 



i=5 



have neither collinear or soft singularities and therefore may be integrated in 4-dimensions. 
In this case the helicity structure is trivially equal to that expected by helicity conservation. 
These differential cross sections are as follows: 



{s, t, u, h) = e 2 f — ^(l - h)C F "\ ~r (t -\- u) 2 {{5Q 2 - s)s 2 - t 2 - u 



da 2 

i2 i „-i„ n2 fr>2 



2s 2 ((5Q 2 + s)s 2 Q 2 + (Q 2 - 3s)tu) 



sA 5 Vp-A 



pX 2 (t 2 + u 2 



-2s 2 (l2(Q 2 + s 2 )stu + s(t + uf - 2s 2 Q A p) 

+4— {2s 2 (3Q 4 ss 2 + 3s(s - s 2 )tu + 2Q 4 (2Q 2 s + tu) - Q 2 s(Q 2 (t + u)+t 
p I \ 

-{Q 2 + s){t + u) 2 tu} 



2 ' u 2 ) 



where 



A = J(t + u) 2 -4Q 2 s 2 , p = Q 2 + s-s 2 . 



(120) 
(121) 



, da 3 ( 
dtdu 
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J S Z7T 
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(123) 



5 Conclusions 



We have presented the parton-level analytical results for the next-to-leading order non- 
singlet virtual corrections to the Drell-Yan differential cross-section. The dependence of 
the differential cross section on the helicity of the initial state partons is shown explicitly 
(the spins of the final state partons are summed) . Although the calculation is very involved 
due to the presence of chiral projectors in the initial state, the result is quite compact and 
has been presented in a form from which cancelation of the infrared, collinear and infrared 
plus collinear singularities is evident. The non-singlet parton level asymmetry defined by, 



da L L 

dt 



da 



++ 



da _j 



dt 



— s- 



dt 



(124) 



is given by order a 2 by, 



dn LL 

dtdu [S,hU) 



d °NS d °o at 



dtdu 



dtdu 



2s 



daf 
dtdu ' 



where da^ s is the non-singlet unpolarized spin-averaged cross section given by 



daV, daV daV da? 
s-r-r- + s-r-r- + s— -f- + s- 



dtdu 



daf 

dtdu dtdu dtdu dtdu 7 



(125) 



(126) 
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where das is the order a s Born term and the remaining cross sections are of order a 2 . The 
leading order parton level result that 



is not true at order a 2 in the MS scheme. In this scheme, Pqq(z) = AP qq (z) in D-dimensions 
and the helicity non-conserving term, daft 0,1 is explicitly needed in the convolution of the 
hard scattering with the NLO evolved structure functions. The evolution DGLAP kernel 
has obviously to be evaluated in this same scheme. 

Our unpolarized cross section, gotten by setting h = 0, agrees with the previous result 
of Ellis, Martinelli and Petronzio jjj in the non-singlet sector. Our calculation is a first 
step toward the extension of the classical 0(a 2 ) Drell-Yan result to the case of longitudinal 
polarized beams. A complementary discussion of our results and of the methods developed 
by us in the analysis of the virtual corrections will be presented elsewhere 
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Appendix A: 2- and 3-particle Phase Space for Polarized Drell 
Yan 

Let's start considering the virtual contributions to the cross section. The relevant phase 
space integral is given by 

(2vr)"-i(2vrV-^ +(g2 " + P2 ~ h - g)(2vr)". (128) 

It is convenient to introduce light-cone variables (q + , q~, q±), with 

q = q+ n + + q~n~ + q± (129) 



and 



and work in the frame 



n± = i=(l,0 ± ,±l) (130) 



pt=P2=\h; ( 131 ) 



to get 

ps 2 = J |^kir /2 - 2 <%il*(2<? + g- - kil - Q 2 )s(s + t + u - Q 2 )n n ~ 3 (132) 

/Qn— 3 
dq + dq-(2q + q- - Q 2 ) n ' 2 - 2 5(s +t + u- Q 2 )^m (133) 
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with 

•n-3 27r n / 2 ~ 1 



fi n ~ 3 = 2 TT siue?- l - z d6i = — (134) 

Jo f = \ 1 r(n/2-l) 1 ; 



using 



d(t,u) 2s 



and after covariantization 



(135) 



(2 *v - Q 2 ) - (Q 2 -*KQ 2 -")-^ 2 = * ( 136 ) 

s s 

PS, _ f^"" ^; /2 "/ (S+ / + ":° 2) (137) 

7 2s V a / r(n/2 - 1) (2vr)"- 2 v ; 

-w.f WT^ a ' if) (^J i(s + t+u ~ Q2)n{Q2 ' u - t] - (m 

Moving to the contributions with 3 particles in the final state, we introduce the invariants 



s = (pi +P2) 2 

sn = (h + k 2 ) 2 S23 = (fa + k 3 ) 2 
h = (pi - h) 2 ui = (p2 - h) 2 

813 = (h + k 3 ) 2 

*2 = (Pi " fc 2 ) 2 t 3 = (pi " fo) 2 

u 2 = (p 2 -A; 2 ) 2 u 3 = (pa - A; 3 ) 2 . (139) 

Only five of them are independent since 

s = Q 2 - (u 1 + u 2 + u 3 ) = Q 2 - (ii + i 2 + *3) = S12 + s 2 3 + si 3 - Q 2 

si 2 = s + u 3 + t 3 

S23 = 8 + Ui + ti - Q 2 

s 13 = s + u 2 + t 2 . (140) 

A discussion of the modifications which appear in the integration over the phase space 
of the 3 final states, due to this ansatz, is presented in the next section. 

The use of the t'Hooft-Veltman regularization || introduces a dependence of the matrix 
elements on the hat-momenta which requires, in part, a modification of the phase space 
integral which appear in the unpolarized case. 

In the case of unpolarized scattering, the singularities are generated by poles in the 
matrix elements which have the form l/t 3 , l/u 3 , l/(t 3 u 3 ) and similar ones, in multiple 
combinations of them. Multiple poles can be reduced to sums of combinations of double 
poles by using simple identities among all the invariants and by the repeated use of partial 
fractioning. This is by now a well established procedure. In our case we encounter new 
terms of the form 1/t 2 and l/u 3 and new matrix elements containing typical factors of 
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the form k 3 , k 2 , and k2 ■ k 3 at the numerator. Let's discuss for a moment these last terms 
containing hat-momenta. It is obvious that by a suitable choice of the parameterizations 
given by the sets 1, 2 ,3 and 4, (defined in the next section) we are able to reduce to 
the ordinary phase space result given by ( |l46j ) all the matrix elements containing scalar 

products of the form pi ■ kj, k\ ■ kj with i = 1,2, j = 1, 2, 3. Therefore it is possible to 
set to zero, after taking the traces, all the products that contain such combinations of hat- 
momenta. Then, the only matrix elements of hat-momenta which are left and which are 

not set to zero are those containing products of the form k a ■ fcb,with a, b = 1, 2. 
The generic 3-particle phase space integral is given by the formula 

po = fstat f d n h d n k 3 d n k 2 
3 ~ 1 J (27T)"- 1 (27T)"" 1 (27T)™- 1 

5(k 2 )5(k 2 )5(k 2 - Q 2 )(2TT) n 5 n ( Pl + P2 -k 3 -k 2 - ki). (141) 

Care must be taken in integrating over k 2 and k 3 when they are identical partons. In this 
case, an extra factor of 1/2 must be inserted so as not to double count the number of 
identical partons. We have indicated this by including a statistical factor f stat that is equal 
to 1/2 for identical partons, otherwise it is equal to one. We use the notation k\ = q to 
characterize the momentum of the virtual photon and we lump together the momenta k 3 
and k 2 as follows 



PS 3 = (2 ^ )2w „3 J d n k 1 d n k 3 d n k 2 d n k 23 5(kl)5(k 2 2 )5(k 2 1 - Q 2 )5 n { Pl + P2 -k 2 - h) 
S(pi +p 2 -h- k 23 )5 n (k 3 + k 2 - k 23 ) 

= (2 ^ )2w _ 3 J d n qd n k 23 5(q 2 - Q 2 )5 n ( Pl + P2 -q- k 23 ) X PS 2 , (142) 



where 

PS 2 = [ d n k 2 5(k 2 )5((k 23 - k 2 ) 2 ) 



{k 2 3 ) n / 2 - 2 ^™/ 2 -3/2 



f d9 1 sin n - 3 e 1 f dd 2 s\u n - i 9 2 . (143) 
Jo Jo 



2 n ~ 2 r[n/2 - 3/2] Jo 
We obtain 

fstat n/2-3/2 r- 

PSs = (2 ; )2 ,-3 2 ,- 2r[n/2 _ 3/2] J ^ - Q 2 )(pi+P2 - q)Y 2 - 2 x I 1>2 (144) 

where 

I x 2 = r dQ x sin n " 3 9 1 f d6 2 sin™" 4 8 2 . (145) 



We now use the light-cone parameterization of q exactly as in the case of the derivation of 
the 2-to-2 cross section, perform the q±_ integration and switch to t, u variables to finally 
get 

(4TT) 2e f stat ( S 



PSs ~ 2WR U^-g^jJ Jl ' 2 ' (146) 
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where we have used the relation s + t + u = S23 + Q 2 and the covariantization 

O + - n 2 Ut ~ Q 2s 23 

2q rq - Q — * . 



(147) 



A.l: Integration of the Real Emission Diagrams 

In order to integrate over the matrix elements, we need to evaluate the various scalar 
products which appear in such matrix elements, in the cm. frame of the pair (1,2). For 
this purpose we define the functions 



It is easy to show that 



\(x,y,z) 
P[x,y,z] 

E[x, y, z] 



2 2 2 

x + y + z — 2xy — 2yz — 2xz 



\ 1 / 2 {x,y,z) 

2^x~ 

x + y - z 

2^ 



(148) 



We get 



\pi\ = P[S23,P1,U1 
\P2\ = P[s23,Pl,h] 

\h\ = \k 2 \ = P[s23,pi,pi] 

\h 

P°l 

P\ = E[s 2 3,pl,tl] 



S23 ~ Ui 

2^i 
S23 - h 



2^3~ 

2 2l _ V^23 



P[S,Q ,S23\ = 

S23 ^V^23 



s 2 3 - U\ 



2^23 

S23 - h 



Pi 



P2 



2^ 



^23 ~ ^1 

2^iii 
S23 - h 
2^ii 

Q 2 - S + S 23 



^2 



2^ 

/S23 



(149) 



(150) 



In the derivation of ( |15C| ) we have used the relation s + 1\ + ui = Q + S23. There are 
four different parameterizations of the integration momenta which we will be using. In the 
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first one, which is suitable for unpolarized scattering one defines (in the cm. frame of the 
pair (1,2)) in general 



^23 



2 

/«23 



1, cos 02 sin 6*i, cos Q\) 



COS (72 sin #i , — cos t>i j 



px = Pi(l,0, ...,0,sin^i,cos^i) 

p 2 = p 2 ) (l,0, ...0,sin^i,cosV>i) 

k\ = (k±,0, ...0, \ki\ sin^, |^i| cos^); 

where the dots denote the remaining n — 2 polar components. 
It is convenient to use the parameterizations 



set 1 



pi =pS(l,0,..., 0,0,1) 

p 2 = p°(l,0,...,-sin<,0,cos<) 

k\ = (ki,0, — \ki\ sin^,0, \ki\ cosip), 



(151) 



(152) 



set 2 



Pl =p°(l,0, ...,sin^",0,cosV") 

p 2 =p°(l,0,..., 0,0,1) 

ki = (ki, 0, \k\ | sin ip', 0, \k\ | cos if)'), 



(153) 



set 3 



Pi = Pi(l> 0, sin tp, 0, cos ip) 
P2 =P2({l,0,...,-smip',0,cos'ip') 
h = (kl 0,..., 0,0,1^1), 



(154) 



where 0, ... refers to n — 5 components identically zero. It is straightforward to obtain the 
relations 



cos^ 
cos^ = 
cosV'' : 



(S23 - h)(s 2 3 ~ Ul) - 2.S23S 



{Q 2 



(«23 - *l)(s23 
-8 + «23)(S23 



Ul) 



2s 23 (Q 2 -ti 



A 1 / 2 (s, Q 2 , S 2 3)(S23 - Ul) 

(s 2 3 - h)(Q 2 - s + s 23 ) - 2s 2 3(<3 2 - h) 



AV2( S ,Q2, S23 ) 



(155) 



We select a specific set depending upon the form of the hat-momenta in the matrix 
elements. We choose the center of mass of the gluon pair (2,3). Following the notation of 
references [p3, 25] we define 
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1^ = r dO t r d0 2 sin 6%- A (a + 6 cos 0i)~ fc 
Jo Jo 

x (A + Bcos 61 + Csin #i cos 9 2 )~ l , (156) 

where a,b, A, B, C, are functions of the external kinematic variables . 

This expression trivially contains collinear singularities if k > 1, I > 1 and a 2 = b 2 
and/or A 2 = B 2 + C 2 . The singularities are traced back to the emission of massless 
gluons from the initial state and in the final states. In this special case, it is convenient to 
rescale the integral and let the angular variables ip, tp', ip" defined above appear. The cases 
a 2 = b 2 and A 2 + B 2 + C 2 , a 2 + b 2 and A 2 = B 2 + C 2 and a 2 + b 2 and A 2 ^ B 2 + C 2 are 
discussed in ref . |25| . It is easy to figure out from the structure of the matrix elements which 
cases require a four dimensional integration and which, instead, have to be evaluated in n 
dimensions. This procedure is standard lore. Notice that only two independent angular 
variables 9\ and 6 2 appear at the time. 

Matrix elements containing more than 2 integration invariants (for instance \/t 2 u 2 u 3 or 
1/^2^3 u 2 ) have to be partial fractioned by using the Mandelstam relations. Just to quote an 
example, we can reduce these ratios in a form suitable for integration in the 6\ 2 variables 
by partial fractioning 



1 1 / u 2 + u 3 



t 2 u 2 u 3 t 2 u 2 u 3 \Q 2 - s - u\ 
where we have used the identity 

u 2 + u 3 



(157) 



Q 2 -s- u r l - (158) 

If we let /[.] denote the corresponding angular integral, we get 

I[l/t 2 u 2 u 3 ] = -— ^ (I[l/t 2 u 2 ] + I[l/t 2 u 3 }) . (159) 

Q A — s — u\ 

The integrals on the rhs of this equation are now in the standard 9\ 9 2 form. 
Defining 

r , n m sin™" 3 6i sin™" 4 9 2 , . 



Jo 



[1 — cos 6iY(l T cos x cos 9\ =p sin \ cos 6 2 sin#i) 



we get 



KK-] 


= C[n, i. 


3W: 


j,n/2 


V n [+,-] 


= C[n, i. 


jW: 


j,n/2 


V n [-,+] 


= Vnl~, 






V n [+,+] 


= C[n, i. 


3W: 


j,n/2 



2 J 

2 + sin|] 
2 J ' 



(161) 



where 



C { n, i,Jl , 2--% r '"/ 2 -, 1 - f "/ 2 7 1 - '] J!'" - \ . (162) 

T[n — 2 — i — j\ T z [n/2 — 1J 
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A. 2: Integration Over Hat-momenta 

Let's now consider the modifications induced by the HVBM prescription in the calculation 
of the real emission diagrams. 

PS 2 = J d n k 2 5(k 2 2 )5((k 23 - k 2 f)t 2 . (163) 

We have set k 2 = (k 2 , k 2 ), with 

^2 = ^(l) cos 63 s'm9 2 sin 0i, cos 2 sin #i,cos 9%) (164) 
being the 4-dimensional part of k 2 . We easily get 

^2 

k 2 = ^j- sin 2 B 3 sin 2 9 l sin 2 6 2 . (165) 
Therefore the usual angular integration measure 

n-2 

d^-V = Y[ sin n ~ 1 - 2 6id9i (166) 
i=i 

is effectively modified to 

3 n-2 

dQ, {n -^ = H sin 11 - 1 OxdOi x J] sin™-'- 2 fyc% (167) 

!=1 Z=4 

The intermediate steps of the evaluation are similar to those in the previous section. An 
analogous treatment can be found in ||, to which we refer for more details. 
In n = 4 — 2e dimensions we get 

Therefore, in the HVBM scheme, additional contributions are generated by the inte- 
gration over the hat-momenta in the 3 parton phase space, which are not present in other 
regularization schemes. These contributions are either finite or of order 1/e and therefore 
modify the subtraction of the collinear terms in the total (real + virtual) cross section. 

Appendix B: Some matrix elements 



The list of the integrals relevant for the real emission diagrams can be found in ref [25]. We 
have independently checked all those integrals which are needed for our calculation and we 
have found agreement with the authors. Except for the integrals requiring an n-dimensional 
integration measure the remaining integrals are quite straightforward. 

In our calculation we also encounter matrix elements of the form t 2 /s\ 2 , t 2 j s\ 2 and 

43/2 
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Using the reference frames discussed in this paper it is possible to reduce the angular 
integration of these matrix elements to integrals of the form 

Y<^ = [\mer 4 r A 5 {1 ~ C ° S ° l)k a TTf. (169) 

n Jo Jo (A + B cos 61 + C sm 0i cos 9 2 ) 1 

We can set n = A since - due to the fact that the photon is massive - the integration is 
finite. For instance we get 

(1,1) _ 2Bb (aB 2 - ABb + aC 2 ) ln X 



B 2 + C 2 (B 2 + C 2 ) 



(1,2) _ 2 (aB 2 - ABb + aC 2 ) Bb\n X 



Y A 



(A* - B 2 - C 2 ) + (B 2 + C 2 ) (b 2 + C 2 ) 
(2,1) _ b (Aa B 3 - 2 A B 2 b + Aa B C 2 + Ab C 2 ) 



2 



: (B 2 + C 2 f 



+ [2a z B — A Aa B b + 2 A B b + Aa B z C z 

-AAaBbC 2 - A 2 b 2 C 2 + B 2 b 2 C 2 + 2 a 2 C A + b 2 C A 
1 



2(B 2 + C 2 ) 5 /2 



lnx 



x A-VB' + C { ' 

and similar expressions for the remaining ones. 

The expression of in the cases where collinear singularities appear, is obtained by 
comparing the angular integral in the imaginary part (the s-channel cut) of the box diagram 
to the imaginary part of the Feynman parameterization of the same diagram [ 24 ] . 
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Figures 




Figure 3: Diagrams which contribute to the process q + q^rf+q + q 
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Figure 4: Diagrams which contribute to the process q + q^j* + G + G 
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Figure 6: Diagrams which contribute to the process q + q — ► 7* + G 
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Figure 7: Diagrams which contribute to the process q + q — > 7* + G 
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